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FINAL REPORT 


October 1, 1969 - December 31, 1974 
Radiation Effects Control - Eyes, Skin 
Contract No. NAS 9-9053 


I . Summary 

A study was undertaken to evaluate the possible adverse 
effects on the lens of the eye and the skin due to exposure to 
proton radiation during manned space flight. 

The following information was considered of prime importance 
in the evaluation: 

1. Relative effects of different proton energies. 

2. Relative effects of different irradiation regimes. 

3. Relative effects of different radiation doses. 

4. Relative latent period. 

Consistent with limitations inherent within available 
facilities and experimental design, every effort was made to 
simulate actual proton irradiation which might be encountered in 
space . 

The New Zealand White rabbit and Beagle dog were chosen as 
experimental models. 1 

Proton energies for exposure included 10, 20, 30, 40 and 50 

MeV. 

Total doses received ranged from 37.5 rad to 2000.0 rad. 
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Irradiation regimes included single acute exposures , daily 
fractionated exposures and weekly fractionated exposures* 

Animals were exposed and then maintained and examined 
periodically until data sufficient to meet the objective were 
obtained. 

No significant skin effects were noted. 

Ocular effects varied with the varying exposure conditions. 
However, in cases where induction of lens abnormalities were 
noted , the abnormalities were minimal in severity. Throughout 
the investigation , no animal exhibited serious sight impairment. 
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II . Facility 

The radiation source used in the research effort was the 88 
inch sector focused cyclotron, located at Texas A&M University. 
The cyclotron has a capability of accelerating protons to 50 MeV, 
deuterons to 55 MeV and alpha particles to 120 MeV. Within the 
physical facility there are four irradiation caves, the main 
cyclotron vault, and associated laboratory areas. The areas of 
interest to this project are indicated in Figure II-l with the 
cyclotron control room indicated by #1, the cyclotron vault by 
#2, the cave used for the irradiations by #3, the low level 
counting room, where all dosimetry information was recorded #4 , 
and the staging area for animal preparation outside of the cave 
area by #5. The star within Cave 3 indicates the position at 
which the beamline was terminated, the transmission chamber 
mounted, and the animal irradiations were performed. 

The beam is transported from the cyclotron in vacuum through 
three inch diameter pipe. To bring the beam into Cave 3, the 
protons were diverted out of the cyclotron along the main beam 
line, magnetically bent through the 155 degree analyzing magnet, 
brought down the far wall of the vault area through a switching 
magnet and into Cave 3. The total distance traveled was in ex- 
cess of 150 feet. At intervals along the beam line sets of quad- 
rapole magnets are located for focusing the charged particle beam. 
In addition, sets of slits to restrict the size of the beam and 
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large carbon blocks, Faraday cups, to totally block the beam are 
also positioned at several locations. 

At the termination of the beam line, a mylar film was used 
as the vacuum isolation barrier. The incident beam thus passed 
through the mylar film, the nitrogen gas and polycarbonate film 
electrodes of the transmission chamber, and several inches of 
air prior to intercepting the animal being irradiated. The de- 
tails of the animal setup at this location will be discussed in 
later sections of this report. 
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Ill . Dosimetry 

Physical Equipment 

The transmission ionization chamber used to monitor the beam 
during irradiation is illustrated in Figures III-l and III-2. 

The chamber is essentially two parallel plate ionization chambers 
sharing a common central high voltage electrode and having a 
guard electrode at each extremity. The plate walls are polycar- 
bonate plastic film 2 microns thick with 280 angstroms of aluminum 
vacuum evaporated on them and the gas medium is dry nitrogen sup- 
plied at 0.1 liters per minute. The aluminized conductivity 
coating is broken in one case by a vertical insulation barrier 
and in the second ionization chamber by a horizontal insulation 
barrier of non-aluminized polycarbonate film, thus dividing the 
two parallel plate chambers into actually four ionization currents 
an upper and lower from one chamber and a left and a right from 
the other chamber. These four signals are then brought out to a 
beam balance monitor. Figure III-3 which essentially took the 
signals and summed them (the four of them) to achieve a total 
signal which was representative of beam intensity. At the same 
time a difference signal between the left and the right halves 
of the one ionization chamber was developed and this was illus- 
trative of the degree to which the beam was uncentered - either 
left or right. An error signal was also developed for the upper 
and lower ionization chamber indicating out-of-balance in that 
direction. These two difference signals were minimized by 
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Figure III-3 
Circuit Diagram s'' 
Transmission Chamber Beam Balance Monitor 


V 
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magnetically focusing the beam while at the same time maintaining 
the proper intensity as indicated by the sum of all four signals. 
Thus, when "balanced" the beam was assumed centered in the beam 
pipe and having equal intensity in each of the four quadrants. 

To verify this assumption, every beam was recorded using polaroid 
film and the uniformity of the image on the exposed film was 
visually interpreted to indicate the degree of uniformity. 

A more quantitative measurement of uniformity was also made 
using lithium fluoride rods in a tissue-equivalent mosaic holder - 
as illustrated in Figure III-4. The rods are 1 mm in diameter x 
6 mm long and while they have limitations in their use for charged 
particle dosimetry, when placed in a common plane perpendicular 
to the proton beam they are excellent for indicating uniformity 
of the beam throughout that plane. 

The absolute calibration of the beam was determined using a 
tissue equivalent* extrapolation chamber. This chamber also used 
the very thin aluminized polycarbonate film as an entrance window 
and is similar in nature to a chamber developed by EG&G and re- 
ported in their technical report S-439-R November 1968, Figure 
III-5. The central electrode is one-half inch in diameter with 
a guard ring about this and an overall chamber diameter of approxi 
mately one and one-half inches. The gas used in the chamber was 
tissue equivalent gas composed of 64.4% methane, 32.4% C0 2 and 
3.2% nitrogen by partial pressure. 

*Shonka A-150 Muscle Equivalent Plastic 



Figure III-4 

LiF Rod Mosaic Holder - !] 
Equivalent Plastic Compos 
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1 - TE PLASTIC ELECTRODE 

2 - ELECTRODE INSULATOR - POLYETHYLENE 
3 - TE PLASTIC GUARD RING 

4 •- HIGH VOLTAGE INSULATOR 

5 - ALUMINUM RING CONNECTOR 

6 - ALUMINUM VARIABLE SPACING RING 



7 - ALUMINUM WINDOW RETAINING RING 

8 - ALUMINUM TUBING CHAMBER STEM 

9 - POLYSTYRENE SIGNAL INSULATOR 

10 - UG-260/U CONNECTOR 

1 1 - HIGH VOLTAGE WIRE 

12 - 0.0005" ALUMINIZED MYLAR 
13- GAS PORTS 


Figure IH-5 

Tissue Equivalent Extrapolation 
Chamber - Cross Sectional View 
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Total particle fluence was determined by a Faraday cup which 
was a disc-shaped piece of graphite located at the point at which 
the animal would be placed and electrically insulated from its 
surroundings . 

A block diagram of the proton spectrometer system is illus- 
trated in Figure III-6. The system is composed of commercially 
available equipment and was used to obtain the proton spectra as 
it penetrated through various depths in the tissue equivalent 
plastic. For the higher energy protons used, the range of the 
proton exceeded the sensitive thickness of a single detector. 

For these measurements a series of three detectors were run in 
parallel into one single preamp and thus into the amplifier sys- 
tem and the rest of the system. This technique maximized the 
problem of the detector noise but minimized cost in that only one 
preamplifier and amplifier were required as opposed to three pre- 
amplifiers, three amplifiers and a summing amplifier, which would 
have been required, had a separate amplification system been 
developed for each detector. The latter technique would have 
minimized the total system electronic noise but for the energies 
used in this project, noise was not a problem. 

The thermoluminescent readout system used for proton beam 
uniformity checks was a commercially available TLD readout system 
and was used in the conventional manner for this project. 



Si (Li) 
Detectors 

Protons I 2 3 

n r~i r~i j£ 
Junction 
Box 



Figure III-6 

Proton Spectrometer System 
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Equipment Characteristics 

Prior to animal irradiations the operational characteristics 
were determined when the equipment was used in detection of pro- 
ton beams. Saturation curves were first developed for both the 
transmission chamber and the tissue equivalent chamber. These 
curves are shown in Figures III-7 and III-8. Based on this data 

^ 7 i 

an upper limit of 8 x 10 amps was chosen for the maximum current i 

I 

which could be produced with the transmission chamber and not ; 

] 

have ion recombination problems. This limit then placed an upper j 

l 

limit on the proton beam intensity and the dose rate which could • 
be achieved. A similar decision was made concerning the extra- 
polation chamber and this limit correlated well with the limit 
that had been decided on for the transmission chamber. For the 
transmission chamber to accurately determine the dose delivered, 
a linear correlation was required with the Faraday cup and the 
T.E. extrapolation chamber; which were used for the absolute 
dosimetry measurements. These correlations are shown in Figures 
II1-9 , III-10, III-ll. When applicable, the current limits im- 
posed by the saturation curves are indicated. As indicated by 
the graphs, there is a very good linear correlation between each 
pair of detectors . 

Basic Dosimetry Theory 

The absolute dose was determined using the fundamental Bragg 
Gray cavity theory. The computation of dose was based on the use 
of equation 1 below: 











Faraday Cup ; Amps 
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Figure III-9 

Transmission Ionization Chamber 
Current vs Faraday Cup Current 
For 30 MeV Protons 


Saturation Limit on TC Current 
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Figure III- 10 

Tissue Equivalent Ionization Chamber 
Current vs Transmission Chamber 
Current for 30 MeV Protons 
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Figure III-ll 

Tissue Equivalent Ionization Chamber 
Current vs Faraday Cup Current for 

30 MeV Protons 
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eq 1 


[m Kg] 


where 

f • 

W 

e 

* 


= Joules/coul for formation of ions = 30.5 


[■ 


= Stopping power ratio (variable with proton energy) 


TE plastic 
TE gas 

4 

[m] = Kg of gas in sensitive volume of ionization chamber 
= (.3346 cm 3 ) (1.138 x 10~ 6 Kg/cm 3 ) = 3.8077 x 10 -7 Kg @ 
0°C and 760 mm Hg 

JtpcJ = Temperature pressure correction to that existing when 
chamber was used. 

[ Q ]= Charge collected from ion chamber. 


The gas used in the ion chamber was a "tissue equivalent" mixture 
composed of 3 2.4% CC ^ , 3.2% Nj , and 64.4% CH^ by partial pressure. 
The basic parameters involved in the computation for the various 
proton energies are listed in Table III-l. 

The electrometer used to measure the charge collected was 
calibrated at the Regional Calibration Laboratory in Houston and 
is traceable to the National Bureau of Standards. The mass of 
gas in the sensitive volume of the chamber was determined using 
the measured volume and an accepted density for tissue equivalent 
gas. The absolute dose as calculated for each of the cyclotron 
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TABLE III-l 

- Energy Dependent Dosimetric 

Parameters 

Proton 

Energy 

MeV 

TE Plastic/TE Gas 
Stopping Power 
Ratio 

W/e (J/C) 

10 

1.0107 

30.5 

20 

1.0087 

30.5 

30 

1.0084 

30.5 

40 

1.0075 

30.5 

50 

1.0079 

30.5 
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runs was then related to a given reading on the transmission ion 
chamber. From the previously demonstrated linearity between the 
two chambers any dose could then be delivered by simply taking 
the proper ratio of the transmission chamber readings to the 
given dose in the standard calibration of the day. 

The Linear Energy Transfer, LET, of the beam was measured 
at finite incremental steps and compared to theoretical expression 
given by equation 2. 


dE 

dx 


4 tt z z e 1 * NZ 

~ nPv 7 



,2 m v 



eq 2 


with the following symbol definitions: 

z = Atomic number of charged particle (1 for proton) 

e = Electron charge 

N = Atom density of absorber 

Z = Average atomic number of absorber atoms 

m ® Electron rest mass 

v = Velocity of charged particle 

I = Average excitation energy of the absorber 

g = v /c = ratio of charged particle velocity to speed of light 
r 

i = Shell correction factor 


This equation has been programmed to allow computation of^either 
a mean range or for generating dE vs dx versus depth curves for 
the various energy protons and the materials of interest. The 
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program is specific in nature in that it does allow incorporation 
of the absorbing effects of the transmission ionization chamber 
and the mylar vacuum isolation film in its analysis of the effec- 
tive energy which reaches the animal. A complete listing of the 
program is included as Appendix 1. 

Beam Spread 

Although the beam had been brought some 100 feet through an 
evacuated pipe and should be at this point and time essentially 
parallel in nature, once brought into air nuclear scattering 
occurred and the beam diverged slightly with distance. This ef- 
fect is illustrated in Figure III-12 for 30 MeV protons. To be 
sure that this problem did not hamper us in our research, the 
absolute calibration of the beam for each run was performed at 
the distance the animals were going to be positioned for their 
irradiation. Typically this distance was approximately 3 inches 
from the transmission chamber. 

Beam Uniformity and Alignment 

Beam uniformity was determined in a "macroscopic" nature 
using both the Beam Balance Monitor and single sheet polaroid 
film. For a detailed measurement the LiF rods in the T.E. mosaic 
holder were used. Initially magnetic focusing was utilized to 
obtain the spots of desired size, approximately one inch in dia- 
meter. The thought behind this being to minimize secondary par- 
ticle production and beam contamination due to scatter with the 
beam pipe or collimators. The results of several measurements 




Figure III-12 

Proton Dose Rate Variation With 
Distance in Air from the Transmission 
Chamber 
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for proton measurements made at two different energies and in 
depth are shown in Tables III-2 and III-3. Inspection of the 
tables quickly indicates that throughout the area exposed to the 
eye there are serious non-uniformities with variations of as much 
as 40% in some areas. Attempts were made to correct this by var- 
ious focusing techniques, but with a magnetically focused proton 
beam, it was not possible to achieve the uniformity we felt was 
required. 

The magnetic focusing of the beam was then abandoned and 
uniformity was achieved by defocusing the beam till it had a 
cross sectional ■'area of four to five inches in diameter. Using 
aluminum collimators the central portion of the beam was selected 
and passed through the transmission chamber. For the beam ob- 
tained in this manner, a series of studies were made to determine 
if the scattered radiation would pose a problem. It was found 
that at the energy levels up to 30 MeV, the scattered radiation 
was minimal. The uniformity of the defocused beam was quite ac- 
ceptable being plus or minus 5% throughout the whole area of one 
inch diameter spot, or one and one-half inch spot as the case 
may be, that was used to irradiate the eye. The same defocused 
condition was used for the skin irradiations. In this case a 
two and a half inch diameter circle was used for irradiation and 
Tables III-4 and III-5 indicate the relative dosage throughout 
that circle and also illustrate the extent of the beam scattering 
and enlargement as it passed through air before it was incident 



TABLE III-2 - Relative dose distribution at dog eye; 
Magnetically Focused Beam 
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20 MeV Protons 

- 28 January 

1971 


Zero Depth 





0.989 

1.18 

0.663 

1.07 

1.10 

0.830 

0.884 

0.881 

1.38 

0.999 

0.036" Depth 





0.656 

1.000 

1.26 

1.28 

1.09 

1.05 

0.570 

1.000 

0.891 

1.20 

0.066" Depth 





0.988 

1.18 

0.961 

0.728 

1.22 

0.541 

0.714 

1.14 . 

1.27 

1.11 

0.101" Depth 





— — 

0.904 

0.932 

1.13 

1.05 

0.867 

0.850 

— 

1.18 

1.02 

0.115" Depth 





1.26 

1.09 

.681 

1.47 

1.41 

. 890 

1.14 

.536 

1.08 

1.11 

0.130" Depth 





0.884 

0.796 

1.09 

1.10 

1.25 

0.963 

1.04 

1.08 

0.894 

1.36 


TABLE III-3 - Relative dose distribution at dog eye 
Magnetically Focused Beam 
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30 MeV Protons 

-28 January 

1971 


Zero Depth 





0.712 

1.281 

1.275 

0.931 

0.854 

0.675 

0.951 

0.881 

0.684 

0.595 

0.032" Depth 





0.673 

0.949 

1.23 

0.782 

1.04 

0.761 

.900 

1.05 

1.09 

.895 

0.107" Depth 





0.893 

1.17 

1.12 

1.14 

0.860 

0.746 

1.05 

0.796 

0.724 

0.842 

0.148" Depth 





0.971 

0.879 

1.22 

1.10 

1.10 

0.890 

1.19 

0.850 

0.770 

0.963 

0.199" Depth 





0.671 

1.12 

0.963 

1.10 

0.898 

0.853 

0.858 

1.09 

0.864 

0.972 

0.250” Depth 





0.674 

1.37 

0.927 

0.901 

1.00 

0.895 

0.745 

1.12 

0.931 

0.867 

0.276" Depth 





0.397 

1.27 

1.38 

0.629 

0.953 

0.715 

0.925 

1.14 

0.658 

0.853 

0.286" Depth 





0.586 

1.08 

0.912 

0.953 

1.14 

0.808 

1.25 

1.00 

0.811 

0.866 

0.295" Depth 





0.688 

0.984 

0.938 

1.11 

0.703 

0.828 

1.11 

0.859 

1.06 

0.750 


0 
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TABLE III-4 - Measured relative dose distribution at dog skin 



30 

MeV Protons - 13 

May 1971 




.033 

.546 

.786 

.381 

.0286 


.168 

.923 

.889 

.832 

.830 

1.06 

.038 


.366 

1.07 

1.08 

1.01 

.137 


.877 

1.07 

1.03 

1.11 

.926 

.984 

.438 


.713 

.877 

.854 

1.11 

.505 


.569 

.996 

.745 

r- 

o 

9 

1— 1 

.863 

1*00 

.459 


.0304 

.615 

.983 

.880 

.0307 



Readings within 

2h" D 

collimator 

area. 




TABLE III-5 - Measured relative dose distribution at dog skin 


30 MeV Protons - 20 May 1971 



.0505 

.551 

.961 

.140 

.035 


.273 

1.30 

1.21 

. 933 

1.03 

.881 

.0548 


.327 

1.07 

1.00 

1.07 

.0665 


1.10 

1.03 

1.13 

.922 

.984 

.813 

1.07 


.172 

1.15 

.968 

.947 

.152 


.172 

.993 

.801 

1.09 

1.07 

.989 

.291 


.0397 

1,06 

1.09 

.903 

.0449 



Readings within 2h" D collimator area. 
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on the animals. The radiation field cut-off is actually closer 
to a two and three quarter to three inch diameter circle. Again, 
over this large a spot the uniformity is seen to be plus or minus 
15% which is clearly much better than the magnetically focused 
one inch diameter spot was capable of achieving. To further 
assist in the alignment of the animals and keep them centered in 
the spot, two items were used. The first was a series of cross 
hairs constructed from sewing thread which were positioned such 
that they were centered in the beam. This was verified with the 
photographic film analysis of the beam uniformity in that the 
cross hairs offered sufficient absorption such that they showed 
in the photograph. The second item was a light located on a re- 
movable holder upstream in the evacuated beam line with respect 
to the transmission chamber. When the light was swung into the 
center of the beam line, a light beam was transmitted through the 
intersection of the clear insulating lines that separated the 
upper and the lower halves and the left and right halves of the 
two parallel plates in the ion chamber. This light beam was then 
positioned on central axis and shown directly on the spot that 
we wished to be irradiated. The coincidence of the light beam 
with the central axis of the radiation beam was verified by film 
and proved to be a very convenient method of always properly 
aligning the eye or the skin in the beam. 
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Hair Absorption Study 

Because of the lack of observable effects on the skin at the 
dosage used, some concern developed over the amount of energy 
which the protons were depositing in the hair of the dogs. To 
evaluate this, hair from four 36 centimeter square areas was re- 
moved from a dog and weighed to give a value for the mass per 
unit area which the proton beam encountered. An average value 
of 10.8 milligrams per square centimeter plus or minus 17% was 
obtained. Assuming a chemical composition for the hair of 12.5% 
water and 87.5% keratin, the energy deposited in such an absorber 
by incident 30 MeV protons was calculated using the computer 
code. A total energy loss of .02 MeV was calculated for 30 MeV 
protons. Based on calculations of the energy loss in passing 
through the transmission chamber, it is logical to expect slightly 
larger energy losses for the 10 and 20 MeV instant beams. How- 
ever, the magnitude of each would still be too small for this 
absorption to significantly affect the range of the' particles in 
tissue. Thus, the lack of skin effect at the doses which were 
used in this study cannot be attributed to the fact that the 
animals were not shaved and thus a large energy absorption occurred 
in the hair. 

Dose Beyond Bragg Peak 

Using the computer program representation of the LET distri- 
bution, the dose at the Bragg peak and beyond was calculated for 
the energies 20 through 50 MeV. The results are tabulated in 
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Table III-6 . In essence the results can be summarized by saying 
the dose drops rapidly with increasing depth to a value of less 
than 1% of the Bragg peak dose and that any unusual artifacts 
observed in biological organs located at depths much in access 
of the location of the Bragg peak are not due to radiation dosage 
at that location. 

Results 

The incremental LET as measured for 10 MeV protons in T.E. 
plastic is illustrated in Figure III-13, A theoretical computa- 
tion using the computer program previously described is included 
as well to indicate the degree with which the two concur. There 
is some concern as to exactly how high the peak (Bragg peak) truly 
is. With the experimental techniques used, we were unable to 
measure doses at thickness increments smaller than l/1000ths of 
an inch. It is therefore possible that the actual peak goes 
slightly higher than that which has been recorded. 

There is also a philosophical concern as to what actually 
is a Rad when one is working with highly ionized charged parti- 
cles. To be defined accurately, the dose to a given mass would 
have to be rather uniformly distributed such that one can indeed 
quote in Joules per kilogram or ergs per gram and have it mean- 
ingful. With the very narrow Bragg peaks observed, the given 
dose rate does not exist for a gram of material but exists only 
for a very small mass of material. The doses quoted at the Bragg 
peak are based on this AE/AX ratio and is normalized to what 
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TABLE III-6 - Dose Beyond Bragg Peak 


Proton 

Energy 

(MeV) 

Depth 

(cm) 

Dose 

(Bragg Peak = 1.00 

20 

0.3510 

1.0000 


0.3808 

0.0233 


0.4111 

0.0015 

30 

0.7506 

1.0000 


0.7650 

0.352 


0.7749 

0.0019 


0.7980 

0.00024 

40 

1.282 

1.0000 


1.292 

0.254 


1.438 

0.0329 


1.5942 

0.0025 

50 

1.86350 

1.0000 


1.9570 

0.5920 


1.9940 

0.1960 


2.0190 

0.0542 


2.106 

0.045 


2.3488 

0.0020 


Relative' Stopping Power, (dE/dX); Rad 


o Experimental Data 
a Program (Theoretical) Data 


A A 


a a o 


6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 

Depth In Tissue Equivalent Plastic; 10” 3 Inch 

Figure III-13 
dE/dx for 10 MeV Protons 
in TE Plastic 
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would have existed for a gram of material. The dose used in all 
computations has been this Bragg peak dose. Because of the dif- 
ficulty in accurately measuring the Bragg peak dose, the entrance 
dose was measured for each calibration and using the determined 
entrance to peak dose ratios, the equivalent Bragg peak dose was 
calculated. It must be emphasized that this method of reporting 
dose was requested by the original contract monitor and was con- 
trary to what the investigators felt to be a more accurate method 
of reporting dose - namely proton fluence incident on the subject 
animal . 

Proton spectra measured by the lithium silicon spectrometer 
system is indicated in Figure III-14 for the 10 MeV protons. The 
only interesting point to note here is that the peaks %o broaden 
as one approaches the lower energies due to energy loss and inter- 
action processes which being statistical in nature tend to spread 
out the range of energies contained within the proton beam. The 
values of energy versus depth are plotted in Figure III-15 and, 
if one extrapolates this back to zero thickness an initial energy 
leaving the transmission chamber of 8.6 MeV is indicated. This 
difference from the nominal 10 MeV that is quoted is due in part 
to energy loss in penetrating through the mylar vacuum barrier, 
the transmission chamber, the gas gaps within the transmission 
chamber and to some extent a deviation of the actual energy ac- 
celerated on the part of the cyclotron from 10 MeV. 



Counts Relative Scale 








Proton Energy ; MeV 
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Similar results are presented in Figures III-16 thru III-18 
with a 20 MeV proto^ spectra. In this case, the energy leaving 
the transmission chamber was calculated to be 19.07 MeV and the 
same general comments given for 10 MeV would hold for 20 MeV as 
well. 

For the 30 MeV data presented in Figures III-19 thru III-21, 
the energy leaving the transmission chamber was determined to be 
28.7 MeV. The 40 MeV results are presented in Figures III-22 
thru III-24. The energy leaving the transmission chamber was 
38.4 MeV and for the 50 MeV data presented in Figures III-25 
thru III-27, the energy leaving the transmission chamber was 47.4 
MeV. In the latter two cases, the 40 and the 50 MeV data, there 
is some indication of scattered lower energy protons as part of 
the original beam prior to any attenuation through tissue equiva- 
lent plastic. This is thought to be due in part to scattering 
off the aluminum walls of the vacuum beam transport system by the 
beam after it was expanded to allow collimation and transmission 
of the central inch and half or two and half inches as the case 
may be for animal exposure. Scattering which occurs up the beam 
pipe some distance from the transmission chamber can occur and 
still pass through the hole in the collimator and thus be inci- 
dent on the animal. This scatter contribution can be minimized 
by magnetically focusing the beam down to a small diameter but 
then one is back to the unacceptable conditions discussed earlier 
of major hot and cold spots within the irradiation area. 
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Figure III-16 

dE/dx for 20 MeV Protons in 
TE Plastic 
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Figure III-18 

Measured Energy vs Depth in TE 
Plastic - 20 MeV Protons 
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Figure III-21 

Measured Energy vs Depth in 
TE Plastic - 30 MeV Protons 
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Figure III-22 

dE/dx for 40 MeV Protons in 
TE Plastic 
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Figure III-23 
40 MeV Proton Spectra vs 
Depth in TE Plastic 
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Figure III-25 
dE/dx for 50 MeV Protons 
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Measured Energy vs Depth in 
TE Plastic - 50 MeV Protons 
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Because of the energy of the protons one detector was not 
capable of fully absorbing all of the energy and it took three 
transmission chambers placed in line to totally absorb the energy. 
There is a small dead layer to the front and a slightly larger 
one to the rear of each detector and energy lost in these dead 
layers is not reflected in pulse put out by the three detectors 
in parallel. Thus, if one finds a particle that has a very high 
LET right in the region of a dead layer, there is some loss in 
the size of the pulse outputed by the detectors for that parti-* 
cular energy or situation. The spectra, therefore, have to be 
interpreted in their total rather than individual details because 
only when interpreted in total will these minor artifacts in the 

i 

data be compensated for and one can get an accurate impression 
as to what is happening as the beam penetrates with depth. There 
is also the problem with three detectors that they have to be 
physically aligned concentrically and parallel to the proton beam. 
That is not to say that the active detection volumes within each 
of the three detectors are completely concurrent. The peripheral 
edges may or may not extend to the same point in all three detec- 
tors, Thus, it is possible for particles to deposit energy in 
j the active volume of one detector but deposit energy in the non- 
active volume in the subsequent detectors. The resultant spectra 
is thus distorted and the peaks broadened beyond what would 
actually exist. 
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Comparison with Other Published Data 

During the course of the experiment a paper was published 
by a group from the Biological Laboratory, University of California, 
Davis, in which proton dE/dx measurements in Lucite were measured 
for 20, 35 and 45 MeV protons in Lucite. The measuring technique 
was to use TLD chips located at various depths in the Lucite and 
from this to infer that dE/dx at that location. The shapes of 
the curves were contrary to what we had measured and in trying to 
determine why the difference between our data and this published 
data, an analysis was performed of both our technique and their 
technique for determination of dE/dx with depth. The complete 
summary of this analysis is included in Appendix 2. It would 
suffice at this time to state that we feel the technique of mea- 
surement using the lithium fluoride chips is quite adequate for 
determination of the uniformity of the beam at a given depth but 
due to the physical size of the chips and the rapid variation of 
dE/dx with position in the phantom, particularly near* the Bragg 
peak, the technique of using the chips is not a valid one and 
leads to erroneous results. In addition, based on cursory exami- 
nation of the published spectra, one would also have to conclude 
that the beam was contaminated with low energy scattered protons 
similar to problems we had experienced on prior occasions and 
had to redesign our focusing and exposure techniques to avoid. 


(? 
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IV. Biology 

Objective 

The objective of this study was to evaluate the possibility 
of damage to the lenticular structures of the eye and damage to 
the skin which might be induced due to exposure to protracted 
low energy proton radiation encountered during manned space 
flight . 

Methods and Materials 

The Texas A&M University Variable Energy Cyclotron was 
utilized as the source of proton radiation. 

The irradiations were designed to evaluate the ocular and 
skin effects of proton radiation with respect to relative effects 
of different proton energies, different proton doses, different 
irradiation regimes, and latent period for induction of tissue 
damage . 

In keeping with the objective, it would have been desirable 
to perform irradiations utilizing a spectrum of proton energies 
similar to that encountered in space. As it was beyond present 
capability to duplicate such a spectrum and, in addition, diffi- 
cult to interpret results obtained from an energy spectrum, 
discrete, monoenergetic protons were indicated. The energies 
selected were 10, 20, 30, 40 and 50 MeV as these energies were 
within the range of energies of interest in space exploration. 

Total doses received ranged from 37.5 rad to 2000.0 rad. 

O 

Selected doses were based on those necessary to induce lens 
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environment. Also considered were its ophthalmological charac- 
teristics lending well to examination, its relative ease of 
handling and restraint for both irradiation and examination pur- 
poses, and its proven susceptibility to development of radiation 
cataracts (reported to be a greater susceptibility than man) . 
Lastly, previous radiation cataract studies utilizing the rabbit 
gave this animal some advantage for accurate evaluation of 
collected data. Females were used to maintain one sex throughout 
the colony thereby avoiding any possibility of breeding during 
periods in which animals were held in group cages. 

Rabbits were housed in individual rabbit cages in a vivarium. 
Each cage was 46 cm x 46 cm x 76 cm with food containers located 
at one end and water nipples installed high enough to require 
standing located at the other end. This arrangement forced some 
degree of mild exercise. 

Cage pans were steam cleaned daily and provided with litter. 
Individual cages were cleaned on a regular basis. Occasionally, 
during transport to the irradiation facility and holding immedi- 
ately prior to exposure, animals were held in group cages. Such 
transport was in an enclosed air-conditioned van designed speci- 
fically for animal transportation. 

Food and water were provided ad libitum . 

Animal holding rooms within the vivarium were equipped with 
an air conditioning system which provided a mixture of fresh and 
recirculated air. Temperature of the rooms was maintained at 
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opacities from exposures to other types of radiation, i.e., 
primarily neutrons. Since these dosages were in the range pre- 
dicted to be encountered on a space journey, they were satisfac- 
tory for the skin studies also. 

Irradiation regimes for both ocular and skin studies included 
total dose administered as a single irradiation, total dose 
divided into five daily equal irradiations, and total dose 
divided into five weekly equal irradiations. 

Again, in keeping with the objective to simulate space 
travel conditions, chronic or low dose rate irradiation would 
have been desirable. Cyclotron produced chronic irradiations 
were not feasible. Thus, chronic irradiation was approximated by 
fractionating the total dose. Acute studies necessitated the 
use of the single irradiation. 

Experimental Animals and Animal Maintenance : Initially the 

experimental model for this investigation was female rabbits of 
the New Zealand White, breed. These were obtained as immature 
1.5 to 2 kg animals from a reliable supplier. Rabbits were 
examined, including ophthalmic examination, upon arrival to in- 
sure they were healthy and free of ocular defects prior to experi- 
mentation. The animals were held for a 30-day quarantine period 
preceding onset of experimentation. During this period each 
animal was tattooed in both ears for identification purposes. 

The rabbit was used because of its size, constitution, 
availability and proven adaptability to a laboratory colony 
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approximately 22° C. In addition, a slight positive pressure 
was held in the rooms to aid in avoiding contamination from the 
outside environment. 

As a preventive measure, rabbits were treated against ear 
mites at 30-day intervals by topical application of an insecticide 
containing rotenone as the active ingredient. 

As the study progressed, AKC registered quality Beagle dogs, 
both male and female, were included as an additional experimental 
model. This was done to provide an animal model with a longer 
life span than the rabbit, this, in that respect, more closely 
correlating results with man. Also considered was the reported 
finding that the dog exhibits a somewhat lower susceptibility to 
radiation induced cataracts than man. Therefore, it was believed 
that the choice of experimental models had bracketed the supposed 
response of man. 

The Beagle was chosen over other canine breeds for much the 
same reasons as the rabbit was chosen, i.e., size, constitution, 
availability, proven adaptability to colony surroundings, relative 
predictable good temperament, and abundance of previously reported 
biological parameters and data. 

The dogs were obtained from reputable commercial kennels as 
16 to 20 week old puppies. Initial routine canine immunizations, 
deworming and identification tattooing had been administered at 
the kennels prior to shipment. As with the rabbits, dogs were 
examined upon arrival to insure health and freedom from pre- 
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existing ocular or skin defects prior to introduction into the 
experimental program. A quarantine period of at least 30 days 
was also provided for the dogs. 

Dogs were housed in outdoor runs 2.4 m x 6 m. Runs were 
constructed on a concrete slab with chain link fencing. Located 
at the rear of every run was a completely enclosed shelter 1.5 m 
x 2.4 m. The entire front and roof of the enclosure could be 
opened to allow ventilation and cooling during warm weather. 

Individual runs were provided with two self-feeders and a 
self -water er. Two feeders were necessary to prevent undue 
fighting between dominant animals and more submissive run mates 
vying for control of a single feeder. Wooden pallets were sup- 
plied in each run to provide dogs an occasional alternative sur- 
face upon which to rest. 

Each run contained five animals of the same sex. The 
arrangement of sexes precluded the possibility of breeding and 
pregnancy. 

Runs were cleaned daily. Food and water were provided to 
the dogs ad^ libitum . 

Immunizations against canine distemper, hepatitis, lepto- 
spirosis and rabies were administered at recommended intervals. 
Deworming drugs and other routine medications were administered 
to the dogs as indicated to maintain good health. 

Irradiation Procedure ; Throughout the irradiation procedures, 
animals were randomly divided into the various groups. Randomizing 
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was accomplished through the use of a table of random numbers 
and the identification numbers tattooed in each animal's ears. 

For the initial exposures, rabbits were divided into three 
major groups of approximately 75 animals each. These major 
groups received exposure to 10 MeV protons, 20 MeV protons and 
30 MeV protons respectively. Each major energy group was further 
divided into three subgroups of approximately 25 animals each. 

One subgroup received a single acute irradiation; one fractionated 
into five weekly equal irradiations; and one into 5 daily equal 
irradiations. Additionally, each subgroup contained five dosage 
subdivisions of five animals each. These dosage subdivisions 
were composed of animals exposed to total doses of 37.5, 75.0, 
150.0, 300.0 and 600.0 rad respectively. All major groups, sub- 
groups and subdivisions received their respective exposures to 
both the eyes and a target area of the skin. An appropriate 
group of control rabbits was maintained. 

During exposure the rabbit was restrained in a modified 
squeeze box leaving the eye and the skin target area (on the mid- 
line just anterior to the point of the shoulders) clear of any 
obstruction to the radiation beam. Irradiations were always done 
to the eye of the rabbit first, then to the skin target of the 
same animal. No eyelid retracting device was necessary to insure 
that the eye remained open during irradiation. The target, eye 
or skin, was "aimed" at a crosshair apparatus placed at the end 
of the beam transporting tube, the center of which was in exact 
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alignment with the center of the radiation beam itself. This 
alignment was determined prior to exposure by photographing the 
radiation beam as it passed through the crosshairs. Initially a 
beam with a nominal diameter of 2 cm was utilized. 

Throughout irradiation the animal was observed visually via 
a closed— circuit television monitor to be certain no excessive 
movement interf erred with the radiation beam striking the target. 

After evaluating early data obtained from the original 
rabbit irradiations, a group of rabbits totalling 58 animals was 
irradiated with single exposure at the 30 MeV energy. This group 
was divided into five subdivisions, each containing ten animals 
and receiving approximately 150.0 rad, 300.0 rad, 600.0 rad, 

900.0 rad, and 1200.0 rad respectively. A sixth group of eight 
animals received approximately 1500.0 rad. Irradiations were done 
to the eye first on the entire group using a 2 cm diameter beam. 
Skin irradiations were then done using a 5 cm diameter beam. 

This was followed by two additional 30 MeV energy groups of 
about 50 rabbits each. One group was irradiated with five daily 
equal fractions; one to five weekly equal fractions. Each was 
subdivided into 37.5, 75.0, 150.0, 300.0, 600.0, 900.0 and 1200.0 
rad dosage subdivisions. Another group of about 50 rabbits was 
irradiated with a single acute irradiation to 20 MeV protons. 

This group was subdivided into dosage subdivisions identical to 
the two 30 MeV groups. 

As trends in data obtained from rabbit irradiations were 
developed, dogs were added to the protocol. 



60 . 


Initially, a group of 100 dogs was used, the entire group 
being irradiated with 30 MeV protons. Fifty dogs received a 
single acute irradiation and 50 received five weekly equal irra- 
diations. Dosage subdivisions for the ocular exposures were 
62.5, 125.0, 250.0, 500.0 and 1000.0 rad, each containing ten 
animals. Dosage subdivisions for skin irradiations were 1000, 
1200, 1400, 1600 and 2000 rad, each containing ten animals. 

Both the single and the fractionated groups contained these total 
dose subdivisions. An appropriate group of control dogs was 
maintained. 

The irradiation procedure for dogs was much the same as that 
followed for the rabbits. Irradiation was accomplished utilizing 
a combination of chemical and physical restraint. Chemical 
restraint consisted of a combination sedation/ anesthesia produced 
by an intramuscular injection of a synthetic narcotic plus pheno- 
thiazine-type tranquilizer 10 to 20 minutes prior to irradiation. 
Physical restraint consisted of a sling apparatus and sand bag 
weights to properly position the dog and prevent movement. Eyelid 
retractors were necessary for ocular exposures in the dogs. 
Bilateral eye irradiations were performed on each dog followed 
by irradiation of two skin patches, one on the forequarter and 
one on the hindquarter. Again, the crosshair aiming apparatus 
and closed-circuit television were used. 

After evaluating results obtained from initial dog irradia- 
tions, two additional groups of dogs were included. These two 
groups received radiation to the eyes only and none to the skin 
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as sufficient data had been gleaned from all previous skin irra- 
diations to permit the formation of a definite conclusion 
regarding proton radiation effects on the skin. Each of these 
two groups was composed of twenty dogs. One group of twenty was 
exposed to a single ocular irradiation of 40 MeV protons, ten 
exposed to a total dose of 125.0 rad and ten 1000.0 rad. The 
remaining group of twenty was exposed to a single ocular irradia 
tion of 50 MeV protons with identical total dose subdivisions. 

These final irradiations concluded all animal exposures for 
the contract. < 

Examination Procedures 

Immediately following irradiation, animals were examined 
weekly. It became apparent that weekly intervals were a greater 
frequency of examination than necessary and thereafter rabbits 
were examined every 30 days while dogs were examined every 60 
days . 

Skin evaluation was largely by objective means, mainly by 
gross observation. More elaborate evaluating techniques such as 
biopsy and microscopic evaluation, skin conductivity tests or 
skin secretion evaluation were never indicated by the results 
noted grossly. 

The ophthalmic examinations were performed in a darkened 
room using both direct and indirect ophthalmoscopes and a slit 
lamp. Eyes were dilated for examination using 1% atropine solu- 
tion. At the time of examination, the examiner had no knowledge 
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of the radiation energy, radiation dose, or prior finding for 
the various animals. No unusual restraints were necessary for 
examination purposes in either rabbits or dogs. 

Lens opacities were scored according to the following system; 

0 - No opacity 

+ - Barely detectable opacity 
++ - Easily detectable, but not a sight impairing 
opacity 

+++ - Sight impairing, but not sufficient to result 
in blindness 

++++ - Animal blind due to lens opacity 
While admittedly this was a subjective method of scoring opacities, 
there appeared to be no satisfactory objective alternative. 

For the purpose of reporting data in a more interpretative 
method, a numerical system of grading based on the four-plus 
scoring was developed. A numerical value representing degree of 
lens opacities for any given group was obtained simply by dividing 
the total number of plus signs assigned to that group of animals 
by the total number of eyes in the group. Thus, the maximum 
possible degree of lens change was a numerical value of 400 indi- 
cating every eye in that group was determined to be four-plus 
(totally blind). If every eye in a group were determined to be a 
two-plus, the numerical value for that group would have been 200. 

Photography was attempted in an endeavor to document changes 
found within the lens. However, for a variety of reasons 
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photographic results proved to be far from satisfactory and this 
method of verifying findings was abandoned. 

Animals were maintained and examined periodically until data 
sufficient to meet the objectives were obtained. In the rabbits 
this period was approximately 450 days post irradiation. In the 
dogs, this period varied from approximately 700 days to about 
1400 days. 

Results 

With regard to skin effects, no animal, rabbit or dog, sub- 
jected to irradiation ever exhibited any truly significant skin 
alteration throughout the entire study. No epilation, no erythema, 
no edema and no ulceration were ever observed. This held true 
even after initial doses and beam diameters were markedly increased 
in later irradiations. 

The complete experimental design for irradiation of the eyes 
is illustrated in Tables 1 thru 3. With respect to ocular effects, 
the initial groups of rabbits were considered to be a pilot 
study to suggest the optimum way to proceed with further study. 

Data obtained from these groups are illustrated in Figures 1 
thru 25. 

Figure 1 shows the control data for all 3 groups. 

Figures 2, 3 and 4 show the results of 30 MeV proton radia- 
tions delivered in single, daily, and weekly exposures, respec- 
tively. It is to be noted that in all instances doses of 600 
and 300 rad are greater than control values. Furthermore, the 



TABLE 1: Exposure Schedule For Initial Rabbit 


TABLE 1A: 


I 

Single Exposure 


3735 75 150 300 600 

rad 


TABLE IB: 


Single Exposure 

I — 1 — I — ! — I 

37% 75 150 300 600 

rad 


10 MeV Protons 


Five Daily Exposures 

I I I ! 1 

21k 75 150 300 600 

rad 


20 MeV Protons. 


Five Daily Exposures 

1 — H — 1 — 1 

37% 75 150 300 600 

rad 


Eye Studies 


1 

Five Weekly Exposures 


3735 75 150 300 600 

rad 


Five Weekly Exposures 

1 H — T — I 

37% 75 150 300 600 

rad 


CT\ 



TABLE 1 : Continued 


TABLE 1C: 


f 

Single Exposure 


37% 75 150 300 600 

rad 


30 MeV Protons 


Five Daily Exposures 

i H — i — i 

37% 75 150 300 600 

rad 


Five Weekly Exposures 

I I I — I — I 

37% 75 150 300 600 

rad 


ui 



TABLE 2: Exposure Schedule For Subsequent Rabbit Eye Studies 


TABLE 2A : 


20 MeV Protons 


37?5 


Single Exposure 

1 — 1 — I — 1 — r 

75 150 300 600 900 


1200 
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TABLE 3A: Exposure Schedule For Dog Eye Studies 

TABLE 3A: 


30 MeV Protons 



62% 125 250 500 1000 62% 125 250 500 1000 

rad rad 

TABLE 3B: 


40 MeV Protons 

. I 

Single Exposure 



125 rad 1000 


TABLE 3C : 


50 MeV Protons 
Single Exposure 



125 rad 1000 



• Control • 20 MeV group 
■ Control - 10 MeV group 
A Controt- 30 MeV group 


Figure 1 


Average degree of lens 
change for Control Rabbits 






Figure 2 
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Figure 3 
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150 rad dose was greater than control values for daily and weekly 
irradiations. 

Figures 5, 6 and 7 show the results of 20 MeV proton radia- 
tions delivered in single, daily, and weekly exposures, respec- 
tively. It is to be noted again that in all instances of doses 
of 600 and 300 rad are greater than control values. Of particu- 
lar interest is the fact that all doses except 150 rad were above 
control values for animals exposed weekly. 

Figures 8, 9, 10, 11 and 12 show the results of 30 MeV pro- 
tons delivered in single, daily, and weekly exposures at a single 
total dose. These data also reflect the greater effects for 
higher doses. 

Figures 13, 14, 15, 16 and 17 show the results of 20 MeV 
protons delivered in single, daily, and weekly exposures at a 
single total dose. It was believed that two things were suggested 
by this data. First, the relative magnitude of effect is greater 
at 20 MeV than at 30 MeV, Second, at least at the lower doses, 
protracted irradiation, i.e., daily and weekly, is more effective 
in causing lens changes. 

Figures 18, 19 and 20 show the results of 10 MeV protons 
delivered in single, daily, and weekly exposures at the three 
highest single total doses. It was believed that there was no 
difference from control values for 10 MeV protons. 

Figures 21, 22, 23, 24 and 25 show the results obtained by 
combining single, daily, and weekly irradiation data for each 
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Figure 7 
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Figure 9 
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Figure 10 
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Figure 13 
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energy and each dose. These graphs more or less summarize the 
data in all graphs previously presented. 

While this initial data would not withstand statistical 
analysis, it was considered that it did provide some definite 
information regarding the ocular effects of proton irradiation. 

First, it was felt that the data indicated no effect from 
10 MeV protons, and that the effects from 20 MeV protons were 
greater in the rabbit than 30 MeV. This was readily explained 
in the rabbit. The Bragg peak occurred at the following approxi- 
mate depths in the eye : 

10 MeV - 0.76 mm 
20 MeV - 3.5 mm 
30 MeV - 7.6 mm 

These depths indicated that the Bragg peak, or peak ionization, 
occurred respectively in the anterior chamber, lens, and vitreous 
humor of the rabbit eye. Thus, a greater effect at 20 MeV would 
be expected. 

Second, it was more difficult from the data to determine an 
absolute difference in the effects of the different exposure 
regimes. Even so, it was felt the data suggested a greater 
effect when the irradiations were protracted. It was also felt 
that weekly exposures may have been more effective than daily, 
particularly at the lower doses. 

Third, there was little doubt that the data suggested a 
greater effect at the higher doses used. 
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Lastly, the data suggested that a latent period was observed. 
While the relationship between dose and latent period could not 
be determined, it was believed that the latent period at the 
higher doses was between 150 and 200 days. 

The above data provided the information required to proceed 
with the study. 

Figures 26 thru 29 illustrate data obtained in subsequent 
rabbit exposures. 

Figure 26 exhibits data obtained from rabbits exposed to a 
single acute irradiation of 20 MeV protons in total doses ranging 
from 37.5 rad to 1200.0 rad. Every dosage subdivision within 
this group exhibited lens abnormalities greater than control 
throughout the post exposure period of study. This difference 
became less with increasing time post exposure presumably because 
of an increasing degree of lens abnormalities in control eyes 
due to normal aging processes. A dose dependent response was 
not noted. 

Figures 27, 28 and 29 illustrate data obtained from rabbits 
exposed to 30 MeV protons in a single irradiation, five daily 
fractions and five weekly fractions respectively. All three 
groups showed a greater degree of lens abnormalities in exposed 
eyes than in control eyes. Again, this difference became less 
with advancing age in the animals. There was no really signifi- 
cant difference in overall degree of lens changes between the 
single and fractionated exposure groups. Also no really apparent 



Figure 26 


* Control 

■ 37.5 red. 

* 75 rad. 

t 150 red. 

* 300 rad. 

* 600 rad. 

0 900 rad. 

2 1200 rad. 


20 MeV Proton c 
Single Exposure 
Rabbits 



Tims Pott irradiation (Daya| 



• Control 

M 37.5 rad. 

♦ 75 rad. 

A 150 rad. 

* 300 rad. 

* 600 rad- 

• 900 rad. 

X 1200 rad. 
0 1500 rad. 


Figure 27 

30 MoV Protons 
Slngio Exposure 
Rabbits 



V£> 

CO 






Figure 29 
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dose dependent response was noted except at the higher total 
doses received ^(600.0 to - 1500,0 rad) invyWhich the trend of higher 
dose, greater response was suggested. 

Figures 30 thru 33 illustrate ocular data obtained from 
exposed dogs . 

Figure 30 summarizes dogs receiving a single acute exposure 

to 30 MeV protons in total doses ranging from 62.5 rad to 1000.0 

$ 

rad. Over the period of study, degree of lens change in the 
radiates was not significantly greater than that in the control 
animals. 

Figure 31 summarizes data obtained from dogs receiving ex- 


posures to 30 MeV proton divided into five equal weekly fractions 
and total doses ranging from 62.5 rad to 1000.0 rad. Only the 
radiates in the highest total dose subdivision (1000.0 rad) ex- 
hibited lens abnormalities significantly greater than control. 
This trend began at approximately 600 days post exposure and 
continued until termination at 1350 days. 


Figure 32 exhibits data obtained from dogs receiving a sin- 
gle exposure to 40 MeV protons. Two total dose subdivisions are 
included, 125.0 rad and 1000.0 rad. Throughout the examination 
period the 1000.0 rad subdivisions showed a slightly greater inci- 
dence of lens abnormalities than control. This greater incidence 
was of so little magnitude that it might be considered insignifi- 
cant.'' With only two exceptions, the 125.0 rad subdivision exhi- 
bited less lens abnormalities than control. ,, 
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Figure 33 summarizes data from dogs receiving a single ex- 
posure to 50 MeV proton in total doses of 125.0 rad and 1000.0 
rad. Lens abnormalities in the 1000.0 rad subdivision were sig- 
nificantly greater than control throughout the post exposure 
examination period. Those animals receiving 125.0 rad did not 
exhibit ocular changes significantly greater than control. 

Throughout the entire study, no animal, rabbit or dog, ever 
exhibited a fully developed cataract resulting in blindness. 

Most lens abnormalities progressed no further than a two-plus 
scoring grade in severity (easily detectable, but not sight im- 
pairing) with only a minimal of three-plus opacities (sight im- 
pairing but not blind) . 

An accidental finding of some interest did occur during the 
course of investigation. Incident to a proposed study unrelated 
to this project electroencephalograms (EEG's) were performed on 
a small number of dogs receiving proton radiation to the eyes. 
Unexpectedly these animals revealed abnormal EEG wave forms . 

With this finding, EEG's were performed on dogs from each dosage 
subdivision, each proton energy group, and each radiation regime. 
Results were read by two qualified examiners, one of which read 
them blind. 

Certain of the dogs exhibited High Voltage, Slow wave 
Activity (HVSA) patterns in the EEG which is interpreted to in- 
dicate neuronal necrosis and is indicative of abnormal function. 
The abnormal EEG changes appeared to be radiation dose and energy 
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related. It was not known whether the observed changes in EEG 
patterns were regressing or progressing. An attempt to further 
evaluate these findings was then undertaken. 

Twenty additional adult Beagle dogs were added to the study. 
Ten had their eyes irradiated with 50 MeV protons to a dosage 
of 1000 rad. Pre-irradiation EEC’s in all test animals were 
interpreted as normal immature patterns as evidenced by mild to 
moderate High Voltage , Slow Activity (HVSA) in all leads. EEG's 
performed at 2 weeks and 2, 5, 7 and 10 months post-irradiation 
revealed no abnormalities. The only change seen in the tracings 
was a gradual reduction in voltage and concomitant increase in 
Hertz. Figures 34 thru 39 illustrate this change in one of the 
irradiated dogs. This is the normal expected change as the EEG 
matures with age of the animal. A normal mature canine EEG is 
shown in Figure 40 for comparison. 

Control dogs exhibited the same EEG patterns as irradiated 
dogs and all were considered normal. Figures 41 thru 46 illus- 
trate the EEG changes seen in one of the control dogs. 



Figure 34 - Pre-irradiation screening EEG (Dog #HHEAS2) under 
pentothal anesthesia characterized by moderate 
High Voltage Slow Activity (HVSA) and interpreted as 
a normal immature pattern. Calibration marks in 
t^is and all subsequent figures: 20 yv and 1 second 
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Figure 35 - EEG recording two weeks post-irradiation (Dog #HHEAS2 ) 











Figure 43 - EEG recording on control dog (Dog #QUIZSL6) two 
months post-irradiation 



Figure 44 EEG recording on control dog (Dog #QIZSL6) five 
months post-irradiation 

| 
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V. Conclusions 

In summary then, we have measured the spectra, the energy 
versus depth and the LET versus depth for the proton energies 
used in the irradiations performed, 10, 20, 30, 40 and 50 MeV. 

From this data we must conclude that the doses which were quoted 
are probably accurate to about 5% based on the constants used. 

The secondary effects of absorption by hair, dose beyond the 
Bragg peak, beam spread with distance, beam uniformity, were all 
considered and either found to be insignificant or a method was 
found in which the effect was corrected for. One concern still 
remains and that is the quotation of dose. The requirement that 
the dose at the Bragg peak be the quoted dose and the dose in 
which all exposures were based, is a bit demanding in that it 
does require one to make a rather "microscopic'' measurement and 
use it in a rather "macroscopic" way. As the dose varies rapidly 
in the peak area one can only assume that the dose measured in 
an infinitesimal thickness indeed exists over a mass as large as 
a gram. Within our ability to measure, which was l/1000th of an 
inch of T.E. plastic, we obtained the curves. This is not to say 
that the curve did not go higher somewhere in between the area we 
were able to measure. The results we have quoted would be in 
error to the extent that the peak height was not properly repre- 
sented by the experimental measurements made. At the higher 
energies it is obvious that there are some problems with scattering. 
On earlier dosimetry runs we had found severe scattering at 40 



122 . 


MeV and for this case the peak entrance ratio was but 2. By 
requiring the operator to focus the beam a bit more and being 
selective as to where the beam was defocused, the beam was cleaned 
up considerably and we were then able to get the data as indicated 
in the Figure III-22. With the 50 MeV beam, a very difficult 
beam to achieve with the TAMVEC cyclotron because it pushes the 
machine to one of its limits of capability , there is still some 
evidence of low energy proton scatter contaminating the beam 
proper. One can only say that the measured spectra and the dE/dx 
curve was that for the beam used to irradiate the animals in the 
experiment. Thus, while it may not be truly representative of a 
pure 50 MeV proton beam, it is representative of the proton beam 
used to irradiate the animals in the experiment. 

As with the uncertain method of quoting "rad dose" with 
protons, one may also question the method of scoring lens abnor- 
malities observed. The types of lesions scored may best be cate- 
gorized as small refractile bodies in the lens. These are known 
to be pre-cataractous in congenital cataracts in some breeds of 
dogs; thus, this was the method chosen. As mentioned above, no 
animals were observed to have developed a true cataract at any 
dosage or irradiation schedule used. 

If our method of scoring lens abnormalities and our method 
of quoting proton radiation dosage are acceptable, the radiation 
dosages now considered acceptable in manned space flight should 
be satisfactory in so far as results obtained in these animal 
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models are concerned. Thus, the primary objective of this study 
has been attained. In addition, no untoward skin reactions were 
observed in any irradiated animal. When considered in the light 
of known skin reactions from other types of radiation, this ob- 
servation is somewhat surprising. 

As with all studies of this type, there are unanswered ques- 
tions as well as factors to be considered in future studies, 
should they be indicated. First, we would strongly recommend 
that the dosage problems mentioned above be resolved. There must 
be a better way. It did not, and still does not, seem appropri- 
ate to use the entry dosage. The integral of the specific ioni- 
zation curve also creates problems for dosage quotations. It was 
suggested that integration of the specific ionization curve over 
the body of the lens be used. This, to be accurate, would require 
measuring the lens thickness and depth in the eye for each animal. 
Even though the rabbits and dogs used were fairly uniform, there 
was still a variation in anatomic size. 

A question still to be answered concerns the dosage required 
to cause proton radiation cataracts. Studies at higher dosage 
levels than those used in this study would be required. Even so, 
this would be essential for further studies. 

It would have been helpful if a standard form of radiation 
could have been used for comparison purposes. In other words, 
determination of an RBE would have permitted us to take advantage 
of a cheaper radiation source as well as data from the published 
literature. 
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No information was obtained which might be applicable should 
radiation modifiers be present. For instance, the high oxygen 
content in the space craft environment might well affect the 
radiation response. 

Lastly, we must conclude that the observed EEG changes in 
the dogs whose eyes were irradiated was a chance observation. 

There are very real and good reasons to doubt this, but in the 
limited study conducted, cause and effect could not be demonstrated. 


I 
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APPENDIX 1 

dE/dx Computation Program Listing 

The program included in the following pages essentially 
calculates the LET versus depth or the mean range for^a given 
absorber, depending on the option chosen, for protons of energies 
between 10 - 50 MeV. All of the constants in the theory involved 
are clearly defined in comment cards in the program. The pro- 
gram has several features which are unique to the experiment at 
Texas A&M in that the transmission chamber and the air gap space 
between the end of the transmission chamber and the animal are 
mocked up in the program. However, this option may or may not 
be included depending on options selected in setting up the pro- 
gram. The energy loss mechanism is represented by equation 2 of 
the main report. 



n n n o n n n .'“i n o n n n n r( n n n n o no n n n o no o o n n o n n n no n n r» n n r> no o n 


1 IJCUBL £ PRECISION AV*"X = N, 7BAR, AT»0F M, P? VASS , £L M A S S, NOPM, P*Jl N»E» 

>*X n 4P T Ni X'UPC v .,.nU3i;. . t gT& t 9=T A ,SHFLL tTPDX ,YDT I N, THICK , R ,ADD,X , Y, VSQ, 

"" *SK,ma , 0 ,00,7 INC.SPl^wC.PUOf ftC,eTA?, PT*6,ftHn,4TMpLC»epF» 
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*A<?0? ,ZA°,0ARS, po , A00? 

2 DOUBLE PRECISION E,S,C, 

3 DIMENSION NAXUt, 4TWTI25), 4TOMSI25), ATNn(25), MS>, 001 5 I , 

! *EXENt2'5 ) ,ATpvpR( 2E) 

4 COMMON 0(4,6), 2INC{4,6) ,AVEXENUO) ,ZBAB < 101 ,ATMDENUO> ,RHO< 101 

5 ‘ IN T FG£P TC, SI, 01 




NOTE — FOP EASE IN UTILIZING TH I S PPOfiP A«, ALL COMMENT CARDS SHOULD 
PE READ, ALL BEAD STATEMENTS AND THFIB EXPLANATIONS HAVE BEEN 
SOUNDED SY ASTERISKS FOR THE IJSEPS CONVFN.I ENCE. 



TH r $ IS A GENERALIZED PROO^A M FDR CALCULATING Thf MEAN RANG# GR 

FOR GENEVA 1 * I NG \ DuD* VERSUS DEPTH CU°VE FOP V AP I A B LF' FN=P GY 

pDOTONS IN DIFFERENT ABSORB FR MATE 0 I 4L S. NORM AL I /A T TON TO 
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5-LFC^ MI AND SET SI = 0 T DJ * 0, AND TC - 0. THE MEAN RANGE 

WILL BE t hc final TOTAL' A*$ORBEP THICKNESS TRAVERSED. THE 

TRANSMISSION CHAMBER MAY *c INCLUDED, BUT THE - RESULTING M£ AN 
M F AN PANGF WILL BE FOR THF ENERGY OF THE PROTON LEAVING THE 
TRANSMISSION AND NOT FOR THE INITIAL INPUT PROTON FNfERGY. 


IF A OF OX VERSUS DFPTH CURVE IS OFSI°EO, SET SI = 1 AND 01 ' ■'i, 

r H£ TRANSMISSION CHAFER CAN BE I NCLUOFD , DEPENDING ON THF 
EXPERIMENTAL SETUP, ACCURATE VALUES OF THE ST ANDAR 0 DEV I A T I ON 
AND MEAN RANGE MUST R? IN^U T ED. A ! 9 CAPS MAY RE INCLUDED, I F 
NORMALIZATION rn F X PE° IMFNTAL DATA IS OE SI R EQ, SET Ml = 1 , N2 - 1, 
AND INPUT ACCURATE VALUES for NOR m , ’ r K A I R 1 # AND TKAIP2, 





NOTE- EXPLANATIONS FOfi DT, Sit MI, TC , Ml, N2 , TKMR1, TKAIP2, 

AND NORM ftR F FOUND BELOW. 


fo 

<T\ 


°* V ASS - PROTON MASS IN CP ft MS 
FLM.4SS * ELer.T = «M l ^ASS IN GRAMS 
C = SPEED O r * LIGHT p. fy/sPC 

TO^TN = thc TOTAL ABSORB 00 THICKNESS IN INCHES. 

F UO FAC - FACTOR USED IN NOR MAL I 7 1 N3 • THF THEO^FTICAI DATA TO 

EXPF^ I m r \ T ft L DA^A. IF NO NORMAL J TATI ON IS to TAKE PLACE, 
FUDFAC WILL AUTOMATICALLY BE T-xEN EQUAL TO 1.0 
4VHXFN = A V F g A GF RXCITATJON FNEPGY (IN FV) OF THE ABSORBER 
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READ* fNINt XBARIN, SIGMA 
*** * * ** * $******* * * ** *£* * ** 


************** ********************* ****************************** 
the a • S = NUMBfR np ARS'i&ApOS OF & SPECIFIED THICKNESS TO BE * 

USED IN THE ORCGR4M. ...*.. 

* 

All) = NUMBER DE 40 KILL THICKNESSCS * 

A( -J) - NUMBER OF 70 MILL THICKNESSES * 

AI3) = NUMBER OF 10 HILL THICKNESSES * 

At 4) = NUMBER Of 5 MILL THICKNESS * 

A(5) = NUMBER OF 1 M ILL THICK NfSSFS 

• * 

THE OQ*S = NUMBER OF OIVlSIONS INTO WHICH C ACH ABSORBER SECTION * 

IS TO BE 01 VIOFO. * 

* 

00(11 =• number OF DIVISIONS in 40 MILL SECT TON* 

001?) = NUMBER OF DIVISIONS IN 20 KILL SECTION* 

— -- 00(B) = NUMBER OF DIVISIONS IN 10 MILL SECTION* 

00(4) = NUMBER Of DIVISIONS IN S MILL SECTION* 

00(5) * NUMBER OF DIVISIONS IN 1 KILL SECTION* 


READ, ( M II , 00(1), 1=1,51 
***** ******************* ** 
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******* ****#-- *******-.* **************************** **********”; 

T X A ! P I = THICK Nr SS OF THE MR GA^"T IN' INCHES) BETWEEN THE 

BEAM TUBE l OR TRANSMISSION CHAMBER r IF USED) AND THE 
ABSORBER HOLOE^o 

OTNCl - NtJM nc P OF OTVIS.ICMS INTO WHICH TKAIP1 IS TO PF DIVIDED. 

TKAI^Z = T HICK V CSS 0^ TrtF air GAP <IN INCHES! SETUP EN THE * , 

' HOLD AND THE f*ETE T COP« IT SHOULD 

MCNTICNEn THAT automatic correction WILL BF HADE FOR 
TH f n EC D F A SE IN THE THICKNESS OF THE GAP DUE TO THE 
placement hf VARIOUS thicknesses of abscpbfp material 
^ -tween THE beam t UPF I op transmission chamber, if 
U S F 0 ) A H 0 r HF DETECTOR. HFNCE. at *0 TIME SHOULD 

- -■ -THE AP GAP THICKNESS ( T K A t R 2 ) BE LESS THAN THE TOTAL’ 

ft BSC!* Ft* THICKNESS USED* 

OtNC? = Nt.’M 9 6° OF DIVISIONS lN r O WHICH TKAIR2 IS TO BE DIVIDED. 

^ryrp AN ACCURATE VALUE FOR TKMR1 MUST BE INPUT ED ONLY^.IF 

^ i\ INITIAL AIR GAP ? S TO *F INCLUDED IN THE PROGRAM. AN 

acc<jra t f value for t ka f r 7 must be imputed only when the 

OfTfCTng is BEING UTILIZED. 


PHAD, TKMR.lt OfXCl, T KA!P2, 0IMC7 

- * ^ ** ** $* * a «*#**& * * * ************** * 


ORIGINAL PAGE IS 
OF POOR QUALITY 



51 

52 

53 

54 

55 

56 

57 
53 
59 

69"“ 

61 

62 

63 

64 

Is 5 
66 
67 
63 

69 

70 
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97 
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95 
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IF 10) CALL ABS n R3(«n 

^ TNT f 4VCX F M ! 1 , T ATMOENfMI) 

CALL OUTPUT t Mf ,«= W j .N t XBA3. IN, ST G M A f TC,OI ? $ T »TKAT81 ff TKAIR?»Nl f N2) 

00 1 t * 1 . « * ===sa ~ — 

1 4V5XENU) = U AV c X5Nl n<‘K60210-19)/U0O-07 

IF ( V T * FO* 10 1 1»-VF**=N<10> = <AVE*£NnO)*l»6O2lD~l9>/U0O-07 
ENIN = CNIN* ( 1 #60210-06) 
xrarcm = xftaRtN*?. 54000 
SIGMA = STG-M4/2. S4DOO 

KKK = ! 

TF (TC.EO.l! CALL 73 ANSI FMT N»ELMASS*C ,PRMA$S I 
IF INl.FO.ll CALL FUDGE! ENT N.TKAfRI , T K A I R2 , 01 NCI ,D! NC2 , NORM, 
*N2,FUDFAC,SL**ASS,PRMASS,C,T07IN,KKK1 • ^ 

PRINT 103 


103 FORMAT! IH1 1 

"'on' ? 1*1.5 

IF (ACT UtO.O.OI GO TO ? 

!F (T.FO.l) T HICK = 0.04 00 00 ' * 

if 11.eo.21 thick = 0. o?oooo . : V . 

IF (I.EQ.3! thick = 0. 01 00 0 0 

IF < I .FQ.4I THICK = 0.00*000 ' • 

JP II.pn.51 THICK = 0.001000 ~ ' - — — 

• NARL =|A( I! 

JSMFS. / 00! I ) 

STFPS7/= THICK/001 I) ‘ 

DELXCAF * STCPS7*?, 54000 " ~ 

00 3 J= 1 »WA e L 

~ enin ' ' — — : 

DO 4 K=1,JA«ES 

P = K ...... - 

Ann = ( steps?*pi 

PRTDIS = 1. 00000 •* • - 

VSC = (C**3.0l*n ,0-f { (PPMASS’»*2.01*(C**4.0) 1 I /! <EM!N ♦ PRMASS * 

*{ C**P. 0 1 1 *“? .0 1 i ' - 

n c .TA = vso/i c**?.ai • 



c^a = beta/ I (1 ,0-&fta**2 1**0. 51 

SHELL = F(M! t CTA1 

OEDX * FI^I, el MASS, VSO, SHELL, 3ET41 ' ' • " “ ' ’ " 

EMIN K EMIN D"LxCM*nEDX • . 

"IF IlFNmi.OD-OT/l. 60210-15). LE.l.ODOOl G n TO 5 

4 CONTINUE 

IF rSI.FO.l) PRTOIS = l.OOOO-GtSIGMA.TOTIN, X0AR1N) ' 

TOTIN = Thtjij * THICK 

SPjnNC = OFOX^PPTOIS/I. 60210-06 

DEnx = O^OX/t. 60210-06 

F = E/1 .6021 0- OS • • • ’ 


EE = EN IN/ 1.60210-06 

P D TNT 100, THICK, TOtin, c, pc , PRTOJS, DFOX, SPIONC 
100 FORMAT ( 51 /) ,720, *TH? THICKNESS OF ThF 45S0R3ER SECTION = 

* *, 016.7,' INCHES* ./ ,T?0,» THE TOTAL APlSPROEP THICKNESS = 

* *,014.7 ,* INCHES*,/ , T 20»’THF PROTON ENERGY ENT gR 1 

*NG THE .\1S0P.PER = ’ ,G14. 7, • Mpy» ,/,T?0, * THE PROTON EWER' 

*GV LEAVING the AESOPPFR SECTION = *,G14.7,* ME V • , / , T 20 . * THF FR 

♦ACT JO* OF ORIGIN At PARTICLES LEFT = * ,G! 4. 7 , / ,T?0, 

**DS0X = * , G1 4,7, ■ «E 

*V/C’-S< , / ,T2o, *THE SPECIFIC TONI 7A T ION IN RAOS/MIN = 

'**,016.71 

"'"IF COt.Eo.ll CALL OFTEC ( EN I N, SPIONC , FUDFAC ,PRTDI St'PRMASSVC* 

♦ELYAS S , TOT IN, TKA I P ?, 0! NO >« KKK 1 
7 CONTINUE 


u> 

o 



100 3 fOM TIMllE 

101 2 CONTI wit* 

102 0RIN T 102 ^ 

103 10? PnOMAT( 61 / W T ?0, • 7HF PROTONS HAVE TRAVERSED ALL ABSORBERS RlfT HAVE 

* N/^fli^OT LOST ALL OF THETR ENERGY* A MEAN RANGE DETERMINATION 
*' ♦ 7 t r 20, Tf> DEC* VERSUS OF 0 TH CURVE IS NOT COMPLETE* THE*,/ f T20, 
«»TO r Al APSOR^ER THICKNESS SHOULD BE INCREASED* M 

1Q4 00 to 20 

105 _ TflTIM s TOTIN' + ADD 

1 06 PR = 1. 7623380*05 ♦ <2.173680*04 * Z BAR ( MI H 

107 F = >NIM/l*60?m-06 * 

108 AD02 ^ <RR / f ft HO ( M I ) *2* 54QOO) ) * E ' 

109 toum s TOT IN + ADD 2 

UO PRINT 101* THICK t F, TOT IN 

til 101 FORMAT (ft< /} *720* *THE ABSORBER ENERGY FEIL ft FLOW l *FV WHILE TPAV 
^FRSING THIS ADSORBERS /.T20* 'THE THICKNESS HE THIS FINAL ABSORBER 

* SFCTjnr * *,G14.7,i INCHES 1 , ST?q** t HF PRCTON ENERGY INT 

*0 This FINAL ABSORBER SFCTtON = •',014*7, 1 MF V* , / ,T70 , •THE TOTAL 

* DISTANCE TRAVELED BY THE PROTON = *,01**7,* INCHES*) 

112 20 PRINT 104 ' ' • 

113 104 FORMAT uni) 

ri* ^ stop '. :*■ ~ 

115 ' END . 1 


116 SUBROUTINE TRans(EMIN*FLMA$S,C ,PRMASS) 

7 11 7 DOUBLE PRECISION AVF ypN , A R , AT Mrjpfg, r , FNI N , op MA S $, E LMA$S , C , SHELL » 

*OFO*, F, $,0100, BETA, ETA, VSO,0*Z INC, EE, STEPSZ 
Tie COMMON 0(4 ,61, ZrNC{4,6) , AVF XEN(lO) VZBARdO) ,ATNDSNI 10) 'iRHOCToY ~ ; 

119 DIMENSION MAXI4J 

c V ' ' 

c '■ ■ 

~ C THIS SUBROUTINE INSERTS THE TRANSMISSION CHAMBER IMMEDIATELY 

C AFTER THE BEAM r^BE AND CA LCUL A T ES THE PROTON ENERGY LOSS 

C THROUGH FT* THE ENERGY 'INTO AND THE ENERGY OUT OF THE CHAMBER ~ 

C ARE OU t PUTED. 

c "' " "" * 

c ‘ 

120 c * ENIN . ■ ; .7 “ 

121 on 1 j = i * 6 __ _ 

122 QU,J) * 1*1 0*06. 

123 1 ZINCUfJ) * ncoooo 

124 DO 2 J~l,6 

125 01 2, J) = P.CH-06 

126 2 7 INC < 2* J I * UODOQ 

127 DC ? J = 1 , 2 

123’ 0 < 3, J ) » 1*250-01 : 

129 3 Z INC t 3 , J ) * IOC* ODOO 

130 CD 4 J = 3 , 4 ’ 

131 o<3, J ) * 1.0000 ^ 

132 “ '4 ' ZIN’Cn, J> ^ 200*00000 \Z 

13.3 013,5) - 0*75000 W 

134 Z ! M C < ?• f ) » 150*00000 *" 

i>5 :<4, n a uoD-oi 

■ 136 Z!\-Ci4,l) ^ ! 00c OODOO 

137 1 1 * 6 

138 MAX<2) = 6 

•!39M5 x< 3> = 5 

~ 160 MAXI 4 ) a 1 *“ "*” * ' " 

14 1 DO 5 ! = 1,4 - 

' 2 MI X I • 



i A 3 kqol = ^ a x ( t 3 

144 on 6 J= 1 ,*<OOL 

145 ?TCP52 = (>MI , J) /7.INC< I , J) )*2. 54000 

146 ' K=/Tvf(l,j) ~ 

147 DO 7 L= 1 ,K 

143 VSO = <C**?. 0 )*( 1 . 0 -( ( ( PP M A$S** 2 . 0 ) * ( C** 4. 01 1 1 /( <EN IN *■ PRMASS * 

*(C**2.0 l)**2.0 I) 

147 BFTA = VSQ/(d**?.0) 

150 F 7 A = BETA/( H , 0 - 8 STA**?. 0 > ** 0 . 5 ) . 

“ 151 ” SHELL = S(MI , £ TA ) ~ ” •* T 

152 DEOX = F(MI, FI. MASS, VSO, SHELL, BETA) 

153 ENIN = FNIN-OE DX*STEPS 7 

154 7 CONTINUE 

155 6 CONT fw»jf 

156 5 CONTINUE 

157 E* F/l . 40710-06 - - - - 

158 FE * EN IN/1. 60710-06 

159 PRINT 100 ,F,FF . , 

160 100 EO°H AT ( 6 ( 7 ) ,’ 30 , 'THE ENFRG Y OF THE PROTON ENTERING THF TftANSMISSI 

•ON CHAM 8 FP = ', 014 . 7 ,' MFV' ,/,T? 0 ,»ThE ENERGY OF THE PROTON LEAVIN 

*G THE T SANSM I S SI ON changer = ',G14.7,» MEV) 

162 FNO 

** WARNING** COMMON PL DCK // HAS A D 1 FEEPEN T LENGTH THAN WAS SPECIFIED IN A PREVIOUS SUBPROGRAM GREATER LENGTH USED. 



1 Z£T 



186 SHELL = Sm,RT ft) 

187 OFDX = F ( v I , EL MA SS ,VSQ t SHELL* B FTA ) 

198 3 ENIN = fNlN - nfnX*5 r EPSZ 

189 IF n.EQol) ENT* - ' ‘ " 

IPO 2 . CGNTTNUF 

191 SPICNC - 0^‘C^pc TOI $*FUDFAC /l* 6021D-06 

19? FF ^ PM l r ’ / 1 * 60710-06 

193 ENIN = F 

_ 194 J E * E / 1 « 602 ID- C6 

199 ‘ v FEE = FFF/U6071D-06 ' “ “ 

196 IF <KKK,FO,fW GO TO 4 

197 PRINT 100, f , PFE, EE t S p IQNC 

198 . 100 FORMAT { 6 ( / ) *T ?0 , * THE PROTON ENERGY^ INTO T HE AIR GAP PROCEED I MG TH 

*E OFTFcTOP = SG14*7,* ME V 0 , / , T?0 /' THE ENERGY OF’ THE PROTON ENTER 

*TNG THE DETECTOR - *,G14.7, } ME V % / , T2 0, * THF ENERGY CF THE PROTON 

★LFAVTNG' THE DETECTOR a *,014,7, * *E v % / , T20 r 5 THE SPECIFIC I ON I Z AT I 

★ON IN MEV/SEC IOP RADS/MIN, IF NORMALIZATION IS P p R FORMED) = «, 

*G 1 4, 7 ) 

199 4 CONTINUE 

200 ' ' RETURN • 

201 _ END 

**WARN!We** COMMON BLOCK // HAS A DIFFERENT LENGTH THAN WAS' S PEC 1 F 1 ED I N'A PREVIOUS SUBPROGRAM GREATER L EMGTH USEO 


202 

203 


204 

C 

c 

c' 

c 

C " 

c 

c 

205 

206 

_ 207 

209 

209 

210 
m 
21 2 

213 ; 

214 

215 

216 
21 7 

213 l 

~ 2 1 9 
220 . 

221 

222 

223 

224 ’ ' " *" 

225 1.01 


SUBROUTINE FUDGE I *NI N # tk & I R 1 , T K A I *2 , D I NCI , 0 1 NC 2 , NOR M , N2 , 

♦FUD^AC, ELMASS, PRMASStC » T OT T N , K KK 1 
DOUBLE PP FC IS ION AVEXEN, 28 A R , A Twr>E N , F , ENI N ♦ PR MA $S , F l w A$ S ,C , SHELL , 
*OEDX,F, S,DLPG, t* , E^A , V SO * 0, 7 I r STE PS 7 , S P I ONC f TK A I PI , D INC1 . 

♦FUOFfif, STEP1 ,EF, T KA[P2 ,0 INC 2, TOXIN, NOPM,o ftTOI c. #EE c 

COMMON 0(4# 6 I , 7INC(4,6 i , avEXENC 10) ,ZRAP HO) »ATmOEN< 10) HO { 10 ) 


THIS SUBOOUTINF WILL INSERT AN ATP GAP QF SPECIFIED T HICKNFSS 
INTO THF PROGRAM AND IS CAPABLE OF CALCULATING A FACTOR WHICH 
WILL NORMALIZE THEORETICAL DATA TO EXP ER I ME NT A L DAT A* 


KKK * 0 

v I ^ B 

STEP 1 = TKA!Pl/r>INC1^2*5 4D0 0 * . 

E * FNIN ■ ““ 

F1JDFAC * 1 o 000 GO 
PP T D I S = !«,0M00 
k = d in r i 

DO 1 1=1, K 

VSQ = (C**?.0^(U0-{MP^ftSS^2.0)MC*«4 # 0)n/f(ENtN ♦ PR MASS * 

*(f>* 2 .on**?«on — 

BETA = VSO/ { C* -2 • 0 ) 

ETA = 3 ETA/( t 1 .0-3ETA**? 1**0. 5 ) 

SHELL =■ SIX I, ETA J ’ 

DEDX - F(MT,BL v A$$,y$G, SHELL, BETA! ■ 

FNIN = FN IN - DPDX*$TE P 1 

I F ( N2, EO.0) GO TO 2 ' 

CALL OF T EC f E N I N# S p IDNC , Ft ro F AC * PR TDI S ? PP M AS S ,C , F L MAS S t 
*T3^lV f T*A IP 2 # DINC2tKKK ) 

FUDFAC = NOPM/SPIONC 
E = P/lo6021D-05 
EE = FN I N< / 1 * 60 210-06 

~P p !N7 101, E , EE, NDPM, FUGFAC 

(5f /) t T?0' ’THF PROTON E N F 0 G Y INTO THF ATP GA D BETWEFN ThE B 
POPT AND ARSOorcr HOLDER = *,G14,7t* MF V • , / , T20 , 1 THE ENERGY OF 






226 

227 

223 

229 

230 

23 r 


232 4 

233 

234 

'♦♦WARNING** 

‘ 235 

236 


237 

233 


* THF PRITON LFAVING THIS AIR GAP » f iG14.7,* *6 V • * / ,720 , *THE EXPFR 
*IM£NTAL VALUE USED IN CALCULATING THE NORMALIZING FACTOR * % 

*G 14. 7 f / T T20* *THE NORMAL T7A I NG FACTOR (FUDFAC) = 1 ? G l*. 7 » / • T20 , 

♦ i-NOPMAt T7 A T I ON HAS BEEN PEREORMEDM ‘ w “‘ 

GO TO 4 

2 CONTINUE ‘ ' ' ' 

£ « E/ 1.6021 0-06 , 

EE = FN IN/ 1/60 210-06 " ’ 

PRINT 100. E, Ec ' - 

TOO FORMAT I 5( / ) »T20» • THE PROTON F NFPGY INTO THE AlP GAP BFTWFFN T HE ~ r : “ 

*8EA* PORT AND ABSORBER HOLDER = *,Gl4*7,» «EV« , / ,T2 0» f THE ENERGY 

♦OF THE PROTON LEAVING THIS AIR GAP = 1 *G14 # 7** M£V* t /* T20t * NO NORM 

♦AH2ATI0N IS PERFORMED. ONLY THF INITIAL AIR GAP HAS BEEN CONSIOE 

*PEO* > % - ■ : 

4 KKK - 1 

RETURN * ‘ " ~ ‘ 

END 

Cn^M^i BLOCK // HAS A DIFFER ENT LENGTH THAN WAS $ PEC IF I EO IN A PREVIOUS SUBPROGRAM GREATER LENGTH USED 

SUBPOUT INF ABSORB I Ml) - "" : ' w < ' ^ ‘ ‘ - . 5 

DOUBLE PRECISION ATP^LC t PHO ? AT WT f ATOMS * ATNH f SUM t AM W * ATHOFN* 

~~"*7 BAR ^ x, Y» A T 0 M FR, E XEN t Z * A VE XfN » OLOG rO* 2 INC 

DIMENSION ATWTI25U ATOMS! 251, ATN0I25), FXENI 25 ) , ATOMFR 1 2 5) 

COMMON 0(4,61 , ZINC (4,6 I , AVE XEN I 1 0) » 7.B AR ( 1 0 ) , ATMOENI 101 f RHOIlO) 


C 

C 

c 

“C 

c 

c 

c 

c 

c 

~c 


This subroutine calcula t fs the average excitation energy 

' i AVEXEN ) * THE AVFRAGF atomic NUMBER i 7 BAR ) 9 AND THE ATOM OWSlTV' 
IATMOFN) FOR A MATERIAL NOT LISTED TN T H£ PROGRAM LIBRARY* 


******£ * *** * *;*: ^^^**#^*5^* A* * *^** + #** ^*#**** <t***^** 

ATPMLC * TflE NUMBER O c ATOMS PfR MOLECULE OF THE ABSORBER* * 

NFLMTS = THF‘ NUMBER OF ELEMENTS IN T HE’ ABSORBER. "V" 


• 

C 

C 

C 

C ■ 
c 

RHO = THF DFSJSITY Of 'HF A B SCIR R E R (G/C«3t. 

ATWT = THE 4.T0MTC WEIGHT OF EACH ELEMENT. ' " ' ” 

aTO^fi = THE NU“3ER OF tTOMS OF F*CH PARTICULAR ELEMENT PER 

molecule of the arsosbep. ’ 

1TN0 =. atomic NUMBER OF FaCH ELEMENT TN THE ABSORBER ' 

♦ 

"■ ♦ 

♦ 

* ' ^ *' 

* 

239 

240 

241 

™c 

c 

3F.10, 1 T PML C • 

rfao, nelmts, PHO 

P540. IATWTIJ), 4T0M$( J) , ATNO<J),J=) .NELMTSi 

* 


■C 



45. 


SUM = 0,0000 

on l J “ 1 * NEl. MT s 

ft M W = SU v ^ATW r <J ) ♦ATOMS! J1 

ATNnFNClOl = I 6. 02 5302 3 1* PHO(IO) ♦ATPMLC / AMW 

"zb ad no) = o.odoo — 

X = O.ODOO 
Y = 0.2 000 

D^ ? J = 1 ,NFLMT^ 

ATO^FR(j) “ ATO’J’Sf JI flT9*LC ‘ 

= 7BARM0) ♦ ATnvf R( J )»ATMO! J) 

IF (ATN0IJ)-13. 0000) 45,45/46 “ 

EXEVIJ) *- ATNi? (J)*(12.0D00 ♦ 7. ODOfl /A TNG! J M 
GO TO 3 


: PET 



255 46 FXENIJJ = AT-NOUl* (9.76000 + 6 R . 8D00 / ( ATNO l J ) **1 « 19 II 

256 3 X = X ♦ ATOMFP (J )*ATNO< J l*DLDG (EXEN< J > I 

257 2 Y = Y + 4T0MFR I J )*ATNOI J ) 

>55'"'”" 7 = X/Y r 

259 AVFXENI 10) = ( ?. 7132fll82fl4) **Z 

260 PP TNT, A VFXEN I 1 0 ) 

261 R £ T UP.N 

262 END 

263~ ' "FUNCTION S(Ml,£TA) 

c 
c 

C THtS IS A SUBPROGRAM TO CALCULATE THE SHELL CORRECTION TFRM 

C 

c 

264 DOUBLE PRECISION AVE XEN, ZB AR, ATMDEN, Z INC ,0 , ET A2 » ET A4 ,ET A6* ET Ar S 

265 COMMON Q{4,6 ) , 7lNC(4,6 I , AVFXFNT 10) ,/BARUO) ,ATMOEN{ 101 ,RH0(10! 

266 FT A 2 = 1.00r>90/(ETA**2l 

267 ? Y A4 = ETA2**2 

265 FT A6 = ETA2**? ' 

269 S = (0. 4?23 7 7r>00*ETA2 ♦ 0. 0 304 043 000* E TA4 -0,00 381 0 600 O’* FT A4 ) * 

*11.00-061* tAVFXENTMI > * * 2 » 0 J ♦' 1 3 . «5 801 9D00*ETA2 - 0.16 67929000"”' 

**5TA4 ♦ 0, 0015795 5000* ET A6 I * < 1. 00-09 I*(AV£XEN(M1 ) ** 3.0 > 

270 RETURN " "" 

271 ENO 

♦♦WARNING** COMMON BLOCK tt HAS A DIFFERENT LENGTH THAN WAS SPECIFIED IN A PREVIOUS SUBPROGRAM GREATER LENGTH USED 

"272 FUNCTION F ( M I , £L M ASS , V SO , SH ELL , BETA ) 7 

C 

C 

C THIS IS A SUBPROGRAM TO CALCULATE DE/OX 

c ........ ... 

c 

273 DOUBLE PR EC ISI ON 0 , 7] NC , AVE XEN , ZBAR , AT MOEN ♦ ELM ASS, V SQ.SHFLL'i " 

*R ET 4 , F , 01 QG 

274 CO«vn N o (4,6-1, ZINC 14, 6) r AVEXC-N (101, ZB AP (10) ,ATMOENI 101 ,RHO( 10 > 

2T5 F = (4. 0*3, 1 41 5^27*1 14. 8 02980- 10 1**4, 0 I *ATMDEN t Mil *7 BAR ( Ml l*(OLOG 

*{ ?.OD | 70*ELMASS*VSO/I AVFXE M ( »I )*( 1 .ODOC-VSQ/ (< 2. 997925010 > **2.0 
*U M-tBFTA**?. 0)-SHFLL/ZfiAR (MI) ) ) /( ELMASS*VSO) 

276 RETURN " 

277 ENO 

♦♦WARNING** COMMON BLOCK // HAS A DIFFERENT LENGTH THAN WAS SPECIFIED IN A PREVIOUS SUBPROGRAM GREATER LENGTH USED 

278 FUNCTION G( SIGMA ,TOTIN,XPARIN» 

C 

C 

C THIS IS A SUBPROGRAM FOP. THF CALCULA T ICN OF STRAGGLING ABOUT 

C THE MEAN PANG? XBAR ASSUMING A GAUSSIAN DISTRIBUTION 

C 
C 

279 DOUBLE PRECISION G, S IGY A ,X9 AR I N t TOTT N,Hl ,H2 ,H3,H4, H5 ,H6 , ARG2.Z AP, 

*OSOP r ,T>ABS ' ' ' 

230 . HI =0.7052307840-01 

281 H2 = 0. 4223201 230-01 

252 H3 =0.927057720-02 

2 s 3 • H4 =• 0. 15701430-03 

2 B4 H5 * 0.27646770-03 

355 " ' M6 * 0.47063 30-04 

f‘,9 AC -,2 » (TOTtN-XBARINl/I 1 OSORT I 2 . 0000 M *S IGMA I 

4.7 IF ( DAP. S( APG2 ) .LF. 1. 00-1 4* Af>G 7 = 0.0000 




288 TFIARG71 1,2,3 

289 ' 1 ARG2 = -ARG2 

290 7A n ^lHl*\RG2 ♦ H2 *&rG2**2 + M3*ARG2**3 ♦ H4*ARC2**4 + H5* 
*ARG?**5 + H6MRG?**6> 

291 IF f 2AP .r,E,50. 0) GO TO 6 

292 0 « 0 . 5 D 0 D + 0.5r>0O*<-l.ODO0+( UOOOO/I 1 1. ODOO+2AP>**16. OH ) 

293 GO Tn 5 

294 2 G = 0.5 

295 GO TO 5 

296 'i CONTINUE 

297 7 APa (HI *ARG2 + H?*ARG2**2 * H3*ARG2**3 ♦ H4*ARG2*+4 + H5 * 
*ARG2**5 * H6*ARG2**6 ) 

298 IF ( ZAP.GF.50. 01 GO TO 7 

299 G *.0.5000 + 0 .50 00* ( 1 . 0000- (I . 0000 /< (1 • QDOO+Z A P) **1 6 . 0 >) ) 

300 GO TO 5 4 

301" 6 G * 0.0 ‘ “ ' . ‘ 

302 GO TO 5 

303 7 G * 1.00000 

304 5 CONTINUE 

305 ' RETURN * " . *" 

306 Fno 




307 SUBROUTINE OUTPUT ( Ml ,ENIN, XBAR I N, SI GMA f TC ,01, S I ,TKA IR1 ,TKAlP2 f 

*N 1 , N 2 J ■ V- 

303 DOUBLE PRECISION FNIN, X9 AR I N, SIGMA ,TKAHU , TKAI 92 

309 INTEGER TC, Or, SI ' " 

31 0 IF (MI.EQ.l) PRINT 20 ' 

311 IF (MI.F0.21 PRINT 21 • ' V"“ ~ ‘ ~ ~ 

312 IF fMI.EQ.3l P’INT 22 

913 IF I M t . EC a 4 } PRINT 23 

314 IF IMI.F0.5I PRINT 24 

‘ “ 315 IF («KEQ.M PRfNT 25 

316 IF (H1.E0.7) P» 1 MT 26 

317 IF <*1.80.91 ‘PRINT 27 

313 IF 1MT.FQ.9) PRINT 28 

319 IF (MI.FQ.m PRINT 29 ~ 

3?0 PRINT 30, ENIN, XHARIN, SIGMA 

32! IF (TC.c 0 # o> POINT 34 

322 IF ITC.EG.l) POINT 35 

3?3 r IF CD!. 50.0) 'PRINT 36 

324 . IF (Ola ?Q. 1 1 PRINT 37 

325 TP I SI. FO.O) PRINT 33 

326 IF I S I, FQ.l) P^INT 39 

327 IF (NUfO.l. ANOa N2.S0. I I PRIAT 3l, T KAIRl , T K A I R 2 

329 TF INI. E0.K4ND.N2.EC.CJ PRINT 32, TKAIRl, TKA1R2 

329 IN1.FO.0.AN0.N2.EO.O1 PRINT 33 

330 20 FORMAT PHI , 101/ ),T? 0 , 1 THF ABS^RflF* MATERIAL UNDER CONS I OF R AT I ON I 

' *S AL U M T MJ* M * I 

331 21 Fhc^at ( I HI , 1 0 C / ) , T?0, 1 The 4BS0R3ER MATERIAL UNDER CONS T DER A T I ON I 

*S P n L YC Ac BOMAT C ) 

332 27 . FORMAT M HI , 1 0 < f K T? 0, * ThF ABSORBFR MATERIAL UNDER CONSIDERATION I 

~ ~*S V! T ^GrM (N?>M 

333 . 23 IlHl, 101/1 ,T20, <THE ARSTP3ER MATERIAL UNDER CONSIDERATION I 

*S M Yl A R 1 \ 

334 24 FORMAT (lH’ f 10(/) ,T?0, 1 THE ABSORBER MAT^iai UNDER CONSIDERATION I 
*S TISSUE FOUIVALFNT PLASTIC 1 ) 

335 25 FORMAT ( 1 H 1 , 10 I/KT 70 , 'THE ABSORBER MATERIAL UNDER CCNS I DERAT ION I 

equivalent gasm * ~ 

3 : f; 26. F 0 R M A T I 1 Hi , 10 l / K T 20, 1 THF ABSORBER MATERIAL UNDER CONSIDERATION I 

*S LUCITE * ) 


1 9£T 



33? 27 fOOMftT ( 1H1, 10 t / » ,T20, 'THE ABSORBER MATERIAL UNDER CONSIDERATION 1 

*S AIR ' ) 

338 _ 23 FORMAT ( 1H1 , 101/ ) , T 20, ' THF ABSORBER MATERIAL UNDER CONS I HERAT [ON I 

" * TISSUE • ) 

339 29 FORMAT < IMJ , 10 ( / I , T 2 0, » THF ABSORBER ’ MA T FP l AL UNDER CONS I HE PAT TON I 

*S NOT IN THF PROGRAM l I0PAR Y' 1 

340 30 FORMAT ( /,t20,*THE INPUT ENERGV IS = ' ,GI4.7,' ME V* , // ,T 20, 

♦ •THF Mr- AN RAMOS FOP THE AD SO n B EB MATERIAL IS = » ,014.7, 

♦' INCHE S',//»T?0, 'THE STANDARD DEVIATION FOR THF ABSORBER MATERIAL 
' * IS .a • *G14. 7, • CM',///) 

341 31 FORMAT t 5( / ) , T 20 * ' THF AIR GAPS ARE BEING INCLUDED* ,/»T20t'TKAlRl = 

* *,G14.7,’ INCHES' 720,' TKAIR2 = • , G 1 4. 7 , • I NCHFS « ,/ , T2Q , * NORMAL 
*IM T ION WILL BE PERFORMED') 

342 32 FORMAT! 5(/) ,T?0, 'THE AIR GAPS A P E BRING I NC LUDFD*, / ,T?0, 

*'TKAI0t = ',014.7,' INC HE'S *,/»T7 0»* T KAIP2 = ',G14.7,» INCHES'/, 

' *T20, • NORMALIZATION WILL MOT BE PFRFOP ^rn' ) 

343 33 FnPMAT t 8 ( / I , t 20 , « NO AIR S4PS ARE I NC LUOED » , /, T20, ' NO NORM AL 12 ATIO 

*M IS PERFORMED') 

344 34 FORMAT { /,T20 f 'THE TP AN SM I SS I n N CHAMBER IS NOT INCLUDED IN THE CA 

*L CIJL ATI ON S ' ) 

345 35 FOP MAT ( / , T 20 , ' THF TRANSMISSION CHAMBER IS INCLUDED IN THF CALCUL 

' ' * AT IONS ' ) ' 

346 36 FORMAT ( /,T20,'THE DETECTOR IS NOT INCLUDED IN THE CALCULATIONS' 

*1 

347 37 FOOMAT ! /,T20,'THF D c TE CTOR IS INCLUDED IN THE CALCULATIONS'! 

348 38 FORMAT t /,T20, ' STRAGGLING EFFECTS ARE NEGLECTED') 

349 39 co°MAT ! / , T 20 , • STR A GGL I NG EFFECTS ARE INCLUDED IN THE CALCULATIO 

♦vs' > . : : 

350 RETURN 

3 51 END * 
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APPENDIX 2 

Comparison with Other Published Results 

The comparison of the proton dosimetry results of the Davis 
group with our results has indicated several areas of difference. 
Initially the mean ranges based on their published data were 
compared with our calculations and significant differences of up 
to 95% were found. For instance , the 20 MeV proton range was 
calculated to be 141 mills of Lucite whereas a published value 
of 275 mills is given. For 35 MeV protons a calculated range of 
389 mills was found and a published range of 373, reasonably well 
in agreement. For the 45 MeV protons a calculated range of 612 
mills was found and an experimental value given in a publication 
as 708 mills - again a significant error. We then questioned 
ourselves whether our program was calculating the ranges in Lucite 
properly. To verify the program accuracy, the values of ranges 
for 22, 34 and 46 MeV protons in Lucite were calculated and com- 
pared to those listed in the Tables of Energy Loss and Ranges of 
heavy charged particles by Barkus and Berger, which was prepared 
under contract for NASA by NBS. We agreed with the published 
values of Barkus and Berger within 7 tenths of a percent, in all 
cases. We therefore came to the conclusion that the ranges in 
Lucite for two of the cases are erroneous. 
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Using the computer program we next determined the entrance 
and the exit dE/dx in lithium fluoride chips positioned in Lucite 
for 20, 35 and 45 MeV protons. Using the thickness of 35 mills - 
essentially one millimeter - for the lithium fluoride crystals, 
which closely approximated the thickness of the crystals used by 
the Davis group; the crystals were centered at different depths 
in Lucite that approximated the positions chosen by the Davis 
group as indicated in their paper . The computer runs were then 
made for one chip at a time. This was done to eliminate the 
possible large energy degradation which would occur if many crys- 
tals were lined up in a single row. The dE/dx of the protons at 
the entrance and the exit of the crystals was computed and the 
results plotted in Figures Al thru A3. As can be seen from the 
20 MeV results, Figure Al, the variation in dE/dx across the 
crystals is relatively small until the mean range is approached 
at which point the variation becomes considerable. As is plainly 
visable, the variation in dE/dx across the crystals centered at 
80 mills is relatively large while the dE/dx variation for the 
crystal centered at 100 mills is so large that it is difficult 
to see how a specific value of dE/dx for that single position 
can be obtained. The same results are observed for the 35, and 
45 MeV data. Figures A2 and A3. In all cases, the energy of the 
protons went to zero before completely traversing the last crystal. 
The size of the lithium fluoride crystals used was so large it 
appears to us that there is no way to generate an accurate dE/dx 
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Figure Al 

dE/dx vs Depth in Lucite for 
LiF Crystals - 20 MeV Protons 
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Figure A2 

dE/dx vs Depth in Lucite for 
LiF Crystals - 35 MeV Protons 
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Figure A3 

dE/dx vs Depth in Eucite for 
LiF Crystals - 45 MeV Protons 
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curve in the Bragg peak region using this technique. The Bragg 
peak measurement discrepancy is further illustrated by comparison 
of the depth dose peak full width at half maximum versus proton 
energy in Table Al. The differences here are obvious and are 
consistent with the results one would expect due to the size 
problem of the lithium fluoride. The net result of this Bragg 
peak measurement discrepancy is reflected in the Bragg peak doses 
quoted for the animals irradiated. With the broadening of the 
peaks, a larger fluence of protons are required to give a stated 
Bragg peak dose. Thus, the total energy deposited per unit of 
dose quoted is greater for the peaks which are broadened due to 
contamination by scattered protons. Different biological results 
per given dose would therefore not be surprising. 
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TABLE Al - dE/dx Full Width at Half Maximum 




TAMU 

DAVIS 

RUSSIAN 

10 

MeV 

. 13 mm 

— 

— 

20 

MeV 

. 57 mm 

2 . 3 mm 

— 

30 

MeV 

* 9 mm 

— 

— 

35 

MeV 

-- 

2,2 min 

— 

40 

MeV 

1.0 mm 

— 

-- 

45 

MeV 

— 

7 mm 

— 

70 

MeV 

— 

— 

6 mm 



